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ABSTRACT 

The numerical kernel approach to difference imaging has been implemented and ap- 
plied to gravitational microlensing events observed by the PLANET collaboration. 
The effect of an error in the source-star coordinates is explored and a new algorithm 
is presented for determining the precise coordinates of the microlens in blended events, 
essential for accurate photometry of difference images. It is shown how the photomet- 
ric reference flux need not be measured directly from the reference image but can be 
obtained from measurements of the difference images combined with knowledge of the 
statistical flux uncertainties. The improved performance of the new algorithm, relative 
to ISIS2, is demonstrated. 

Key words: methods:statistical - techniques: image processing - techniques: photo- 
metric 



1 INTRODUCTION 



Over the last 15 years, gravitational microlensing l|Einsteinl 
Il936h has been observed routinely and used in the study 
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of da rk baryonic matter (lAlcock et al.l 19931: Aubourg et al.l 
19931) and stellar atmospheres ("Albro w et al.ll 19991 . l2001al lb[ 



Fields et al.ll2003l : Ic assan ct al. 2004), and in t he search for 
extrasolar planets lAlbrow et al.li 2000. 20o3 iGaudi et al 
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Bond et al 
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2004 
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Udalski et al.. .2005i : iBeaulieu et al 
Gaudi et al.[|2008l ). The PLANET 



collaboration ( Albrow et al. 19981 ) operates a number of 1 
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m class telescopes distributed around the Southern Hemi- 
sphere and performs round-the-clock CCD photometry of 
microlensing events that have been discovered and alerted 
in real time by the OGLE (Udalski et al. 1994; Udalski.2003 ) 
and MOA jSond et al.ll2002l ) microlensing surveys. 

In this paper we discuss recent advances in the 
PLANET difference imaging reduction pipeline, focussing 
on several subtleties inherent in the reduction of blended 
microlensing events. In what follows we use the convention 



D 



R(g)K -T 



(1) 



for difference image D, reference image R, target image T 
and convolution kernel K. That is, a difference image is de- 
fined as the convolved reference minus the target, normalised 
so that it is on the effective exposure scale of the reference. 

We illustrate the methods using a sample dataset of im- 
ages of microlensing event OGLE 2008-BLG-229 that were 
taken using the Elizabeth 1.0-m telescope at the South 
African Astronomical Observatory during PLANET opera- 
tions in 2008. The microlensing event was alerted by OGLE 
on 2008 May 3 and initially was predicted to have low 
magnification. Subsequent observations revealed a blended 
moderate magnification event, peaking with magnification 
Ao = 7.24 on 2008 July 18. OGLE data and parameters for 
the event can be obtained from the OGLE Early Warning 
System website 0. The SAAO observations consist of 84 im- 
ages, spanning the time from 7 days before maximum until 
23 days after maximum. 



2 THE REDUCTION PIPELINE 
2.1 Introduction 

In our first years of operation, PLANET photometry was 
performed both in real-time at the telescopes and offline 
using the DoPHOT PSF-fitting code fec hechter et al.lll993l ') 
under a reduction pipeline written mainly by JPB. Followin g 
the development of the ISIS code l|Alard fc LuptonI [l998l ) . 
we adopted the difference imaging method for obtaining our 
best photometry offline, while still employing the DoPHOT 
pipeline at the telescope sites. In the 2006 season, we began 
using a difference imaging photometric pipeline at the tele- 
scope sites and for our final offline photometry. The offline 
version, known as pySIS2, w as deve l oped by MDA, and is 
based on the ISIS 2 code of I AlardI (|200(]| ). An adaptation 
of pySIS by CC, known as WISIS, is used for most of the 
real-time at-telescope reductions, while the pySIS2 code is 
used at the Perth Observatory. The pipelines allow single 
images to be reduced immediately after observation using 
an existing reference template. As better quality images are 
acquired, the reference template can be updated and the 
previously observed images rereduced. 



2.2 Image registration 

The ISIS code requires that all images be fully registered 
to an astrometric reference. Bright stars are located on all 
frames and cross-correlation of their positions followed by 
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an iterative rejection scheme is used to define an astrometric 
transformation for each target image. 

An innovation introduced to our offline pipeline in 2007 
was the removal of the requirement to fully register images. 
Instead, we register the images only by a shift in X and Y 
to the nearest pixel, thus avoiding the need for interpolation 
and resampling. Resampling is generally undesirable since 
it introduces correlations between adjacent pixels, meaning 
that their ffux uncertainties are no longer described by Pois- 
son statistics. In the case of images that are close to or below 
critical spatial sampling, resampling introduces an artifact 
where stellar PSF's are not constant or slowly-varying across 
an image, but depend on the subpixel location of their cen- 
troids. Such images usually do not subtract cleanly. Integer- 
pixel registration was handled in our modified version of ISIS 
by offsetting the kernel centroid by the subpixel registration 
residual. We note that this approach is somewhat less flexi- 
ble than the standard ISIS code, in that it cannot work with 
sets of images with rotations relative to each other. 

2.3 Difference imaging with a numerical kernel 

In 2008, we have developed a new version of the code, py- 
SIS3, that is no longer based on ISIS image subtraction. In- 
stead, for the di fference-imaging step, we have implemented 
the algorithm of iBramic hl (|2008l ). In this method, the kernel 
is represented as a numerical pixel array, rather than the 
decomposition of Gaussians multiplied by polynomials used 
in ISIS. The numerical kernel is able to accommodate im- 
ages with irregular PSFs, for instance trailed images, that 
ISIS cannot cope with. An implicit feature of the method 
is that complete registration is not required and the kernel 
naturally incorporates subpixel offsets. Image re gistration 
in pySISS is hence restricted to integer pixel shifts. BramichI 
(2008) shows examples of how the new algorithm outper- 
forms ISIS. DMB's code has been used successfully to dis- 
cover new variable stars in the globular cluster NGC 6366 
I Arellano Ferro et al.ll2008l '). 

Our implementation has been used for the analysis of 
several microlensing events, a ppearing in forthcom ing pa- 
pers on MOA 2007-BLG-197 (|Cassan et al.ll2009l ). OGLE 
2004-BLG-482 (IZub et al.ll2009l ). and OGLE 20 07-BLG-472 
jKains et al.l |2008| ). and for a transit search l| Miller et al.l 
12009'). There are several speciflc details of our implementa- 
tion that we note here. 

First, the algorithm is generally more computation- 
ally intensive than ISIS, and the computation time scales 
strongly with the number of pixels required for the kernel 
array. For microlensing events, we generally reduce only a 
subsection of the images, typically 250 x 250 pixels centred 
on the microlens. Our plate scales are typically around 0.3 
arcsec/pixel. 

Second, results depend on the size of the pixel array cho- 
sen for the kernel. The kernel needs to be large enough to en- 
capsulate the transformation between reference and target, 
but not so large that it introduces noise into the convolu- 
tion. We have found the best results by employing a circular 
kernel with a radius (in pixels) given by 

-Rkcrncl = min (7, 4 (FWHMtargot - FWHMrefcronco)) (2) 

where FWHMtargct and FWHMrofcrcncc are the full-width at 
half-maxima (in pixel units) of the microlens on the target 
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and reference images. Regions of the kernel that are located 
more that 7 pixels from its centre are represented by 3x3 
binned pixels in order to reduce noise. The values for the 
binned kernel pixels are computed from the equations in 
iBramichI l|2008h but using a 3 x 3 boxcar-smoothed version 
of the reference image. 

Third, to prevent saturated stars that have irregular 
PSFs from entering the kernel determination, we mask a 
circular area of radius 15 pixels around all pixels that are 
saturated in either the reference or target image as well as 
masking the microlens itself. The default behavior is to use 
all t he remaining im age pixels to determine the kernel using 
the iBramiclil ()2008 ') algorithm. In cases where images are 
contaminated by artifacts, such as diffraction spikes, that are 
not easily masked, we have found that using 'stamps' around 
bright unsaturated stars rather than the entire unmasked 
image renders a kernel that is less prone to contamination. 

Fourth, the photometric scaling factor. 



(3) 



where K is the convolution kernel, represents the relative 
difference in effective exposure time between the reference 
and target images, i.e. it accounts for differences in both ex- 
posure time and atmospheric transparency. For each target, 
if s is significantly different from the ratio of true exposure 
times, this usually indicates a poorly subtracted target or 
one affected by cloud. 

Fifth, for images that have poor spatial sampling - ei- 
ther close to or even below critical sampling - we use the 
following technique. The registered images are oversampled 
by a factor o f two in each di rection using cubic O-MOMS 
interpolation (|B1u et al.ll2001^ ■ This type of resampling does 
not transfer flux across original-pixel boundaries. A stack of 
typically 10 of the best of these images are then mapped onto 
the best-seeing image and combined to make a reference im- 
age. Since our individual images usually have random sub- 
pixel dithers, this process generally results in an oversam- 
pled reference so long as the initial undersampling is not 
too severe. This approach is similar to that employed in R. 
Gilliland's code f or difference-imaging of undersampled HS T 
WFPC2 images l|Gilliland et aLlbOOd : lAlbrow et al.ll2001dh . 
We note that this approach is still under development and 
is not used for the sample data set of images for OGLE 
2008-BLG-229 in this paper, which are not undersampled. 



2.4 Photometry 

To extract photometric measurements from our difference 
images, we first use the Bphot program from ISIS to com- 
pute the PSF of the reference image. This PSF is then con- 
volved with the previously-computed kernel to produce a 
PSF for each target image. The PSF is then normalised 
and resampled at the subpixel lens coordinates using cu- 
bic O-MOMS interpolation. Any residual background is re- 
moved from the difference image using a low-order polyno- 
mial model. Finally, the PSF is fitted to the difference image 
using optimal extraction, i.e each pixel weighted by the in- 
verse of its flux variance. 
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Figure 1. Flux determined at offset coordinates relative to the 
flux at the correct coordinates as a function of coordinate offset, 
AR in units of PSF full- width at half-maximum (FWHM). SoUd 
line for unweighted PSF-fitting, dashed line for optimal PSF- 
fitting with zero background, both under the assumption of a 
Gaussian PSF. The three distinct groups of data points are com- 
puted from the SAAO difference images for OGLE 2008-BLG-229, 
with the PSF shifted in x by 0.2 (green), 0.5 (red) and 1.0 pix- 
els (blue) (0.06, 0.16, 0.31 arcsec respectively) from their correct 
value. 



2.4-1 Flux errors from imprecise coordinates 

All of the microlensing events towards the Galactic Bulge are 
blended to some extent in regular ground-based imaging, i.e. 
the PSF of the microlensed source star overlaps with nearby 
stars. One implication of this is that the coordinates of the 
true source star are often displaced by perhaps several tenths 
of an arcsec from the centroid of the PSF. 

When performing photometry on a set of difference im- 
ages by PSF-fitting, small errors in target location lead to 
systematic underestimates of the fiux. In Figure [T] we show 
how the flux measurement depends on coordinate displace- 
ment, AR (measured in units of FWHM), under a Gaussian 
PSF assumption. The two limiting cases, shown as lines, are 
(a) optimal PSF fitting in the zero-background limit and (b) 
unweighted PSF fitting, essentially the background-limited 
case. In both limits, the fiux error scales as the square of the 
ratio of the coordinate error to the FWHM of the PSF. This 
means that a flux error due to an incorrectly-positioned PSF 
is more serious for images with better seeing and that such a 
coordinate error introduces scatter into a lightcurve derived 
from a set of variable-seeing (or even variable-background) 
images. 

We can demonstrate this effect using our test case of 
SAAO OGLE 2008-BLG-229. We have created difference 
images using as a reference template a single good-seeing 
image, number 18, acquired close to the peak of the mi- 
crolensing event. We have performed optimal photometry 
on the set of difference images at the correct coordinates 
(i.e. zero-offset) and again with the coordinates shifted in X 
by 0.2, 0.5 and 1.0 pixels (0.06, 0.16, 0.31 arcsec) from their 
correct position. The three displayed groups of data points 
in Figure [T] show the flux measured at the offset coordinates 
relative to the flux at the correct position for the three sets of 
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Figure 2. Sample images 78 (left; good seeing) and 83 (right; 
poor seeing) from the SAAO observations of OGLE 2008-BLG- 
229. Top row: direct images registered to the nearest pixel, lin- 
ear greyscale encompassing 95% of pixel values. Second row; dif- 
ference images, linear greyscale range -400 to -1-400. Third row: 
difference images after subtraction of PSF fitted at correct coor- 
dinates. Bottom row: difference images after subtraction of PSF 
fitted at coordinates shifted by +1 pixel in X. 

displaced-coordinate measurements. The figure shows that, 
as predicted, the measured difference flux decreases and its 
dispersion increases with coordinate offset. 

It is important to note that, for difference-image pho- 
tometry, the above effect applies to difference fiuxes, AFi. 
The total fiux is given by Fi — Fq — AFi, where Fq is the 
flux on the reference image, R in Eqn (1). This means that 
the place in the lightcurve where coordinate errors are man- 
ifested most strongly depends on the choice of reference im- 
age. The flux error is largest for regions of the lightcurve 
where the magnification is most difi^erent from that of the 
reference image. Consequently, such errors can be minimized 
for a given part of a lightcurve (for instance some part of a 
lightcurve suspected to display an anomaly) by choosing a 
reference image where the source star has a similar magni- 
fication. 



2.4- S Precise coordinates for blended events 

For high-magnification events, the true source location may 
be discerned from images taken near peak magnification. 



when the flux from the true microlensed source star dom- 
inates that from nearby blended stars. For data sets com- 
prised of images that have precise registration, a sum of the 
absolute values of the difference images can be used success- 
fully to refine the coordinates from their initial estimate, 
even for relatively low-magnification events. This method 
was used in pySIS2. 

In our current circumstance, we have sets of images 
that are registered only to the nearest pixel. The method 
of stacking difference images could, in principle, be applied 
to this data, provided each of the difference images is first 
shifted by the subpixel registration residual. However, such 
a shift requires accurate knowledge of the subpixel residuals 
and involves interpolation and resampling, a process that 
is inaccurate for images with near- or below-critical spatial 
sampling. 

A better way, that retains the original sampling, is to 
use the residuals from PSF-fits to the difference images. 
In appendix!^ we introduce a new algorithm to refine the 
source-star coordinates by minimising these residuals over 
all images. In our photometric code, the algorithm operates 
as an integral part of the measurement process. 

In Figure [2] we show direct, difference and residual im- 
ages for two sample observations, numbered 78 (good seeing) 
and 83 (poor seeing) in the SAAO data set for OGLE 2008- 
BLG-229. Both images were taken during the final days of 
data, when the source was at a magnification, A ~ 2.5. The 
reference image was again image number 18, taken near peak 
magnification, A ~ 7. 

Our coordinate algorithm resulted in a change of 0.54 
pixels (0.17 arcsec) in the location of the target star rela- 
tive to the position found from our best-seeing image (which 
was adopted as the astrometric reference). Using the correct 
coordinates, the residual images (difference images after sub- 
traction of the fltted PSF) are very clean. 

The effect of a coordinate offset of 1 pixel (0.3 arcsec) 
is to produce a residual image with large positive and neg- 
ative flux features remaining. The effect of such an offset 
on the lightcurve can be tested. We compare difference-flux 
lightcurves in a blend-free way by mapping them to a point- 
source point- mass lens model as follows. At time ti, the un- 
blended flux from the source star is given by 

F, = A.Fbase = Fo - AF, (4) 

where Ai is the magnification that we constrain to be de- 
fined by the OGLE geometric parameters for the event (uo 
= 0.139, tE = 53.994 d, to = JD2454665.780), Fo is the un- 
blended source fiux on the reference image and -Fbasc is the 
unblended baseline source flux. For each lightcurve, we solve 
for Fo and Fbasc by minimizing 

2 (^^Fbaso - Fo + AF,f . . 

X (5) 

i 

Lightcurves for the whole data set are shown in Figure O 
for the correct coordinates and for those with a 1-pixel offset. 
An increased scatter is visible in the lightcurve correspond- 
ing to the offset coordinates. In the same flgure, we also show 
the best lightcurves we have derived using a template cre- 
ated from the 10 best-seeing images using pySIS3 and ISIS2. 
Figure |4] shows the corresponding residuals. The effect of a 
coordinate error can be seen through comparison of the up- 
per two panels, particularly during the last 10 days data 
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Figure 3. Lightcurves for OGLE 2008-BLG-229: (a) using sin- 
gle template image from near peak, (b) as for (a) but with a 1 
pixel coordinate offset, (c) combination of 10 images for reference 
template, (d) best lightcurve obtained using ISIS2. 



points when the magnification is most diflerent from that of 
the reference image. The superior performance of the new 
numerical-kernel algorithm can be discerned through com- 
parison of the two lower panels, where the ISIS2 lightcurve 
residuals (panel (d)) have a 53% greater RMS scatter than 
the pySIS3 residuals (panel (c)). We note that there is an ap- 
parent systematic residual in all the displayed lightcurves, 
where the earliest data points lie below the magnification 
curve. This is likely due to the fact that we have constrained, 
rather than fitted, the underlying geometric model. 



2.4-3 Reference flux 

The output of our difference-image photometry is the dif- 
ference in flux, AFi, between each target image, i, and a 
photometric reference image. The reference image may be 
a single observation, as in the preceding sections, but more 
commonly is created from a combination of images with the 
best seeing and lowest sky background. In order to interpret 
our observations we require the flux, Fi = Fq — AFi, where 
Fq is the flux on the reference image. In a model-dependent 
sense, a deblended Fq can be derived as a fitting parameter 
as done above. More usually, particularly during the time 
when observations are being acquired, Fq is measured from 
a PSF-fit at the lens coordinates on the reference image. Of- 
ten this estimation is in error due to the crowded nature of 
the Galactic Bulge fields in which we observe. 

An incorrect value for Fq can, to some extent, be com- 
pensated for in the microlens blend fraction. 
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Figure 4. Lightcurve residuals for OGLE 2008-BLG-229 corre- 
sponding to Fig[3]after subtraction of a best-fit PSPL model with 
the geometric parameters set to be those found from OGLE. Er- 
rorbars are suppressed for clarity. 



/bl 
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(6) 



Fs + Fb 

where for magnification A, the baseline flux {Fs + Fb) in- 
creases to 



F = FsA + Fb = {Fs + Fb) {fbiA + (1 - /bi)) 



(7) 



and the blend fraction is derived from lightcurve fitting. 
However, blending parameters so derived are inconsistent 
between different datasets for the same event, may take on 
non-physical values, and certainly no-longer have the correct 
physical interpretation as the fraction of light contributing 
to the stellar PSF at baseline from the microlensing source 
star. 

A more- successful approach that we have developed is 
to choose Fo so that the photometric uncertainties in the 
set of AFi measurements are consistent with Poisson noise 
for fluxes Fi = Fq — AFi ■ The algorithm for this determina- 
tion, referred to as the Poisson reference flux is detailed in 
Appendix[B] The method uses information from all suitable 
images and the variance in the Poisson reference flux scales 
roughly with the inverse of the number of images. This vari- 
ance is generally smaller than the variance for a direct flux 
measurement from the reference image, which scales with 
the inverse of the number of individual images incorporated 
into the reference. For our sample lightcurve, the reference 
flux is measured directly as Fo = 101, 820 ± 90 ADU (likely 
to include extra blended light), while the poisson method 
yields Fo = 93, 570 ± 40 ADU. 



3 SUMMARY 

Difference imaging has proved to be a powerful technique 
in the measurement of gravitational microlensing flux varia- 
tions and for variable stars and transiting ex trasolar p l anets . 
The numerical kernel method introduced bv lBrarrdc hi (I2008D 
represents a sign i ficant advance over the analytic kernel of 
lAlard fc LuptonI l| 19981 ). 

In this paper we have shown that, for blended microlens- 
ing events, precise determination of the coordinates of the 
microlens is necessary to obtain accurate photometry. We 
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have presented a new algorithm, based on photometric resid- 
uals, to measure such coordinates to high precision. 

Additionally, we have introduced a new method to mea- 
sure the reference flux in a manner that does not depend 
on an (often inaccurate) analysis of the photometric refer- 
ence image. The Poisson reference flux method produces a 
more accurate and precise determination of the unmagnified 
source flux than a direct measurement from the reference 
image. 

A comparison has been made between photometry from 
our new code and photometry using ISIS for SAAO images 
of microlensing event OGLE 2008-BLG-229. The new code 
produces measurements that display significantly less scatter 
about a point-source point-mass-lens lightcurve based on the 
OGLE-determined geometric parameters for the event. 
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APPENDIX A: REFINING THE LENS 
POSITION 

Al target position from a single difference image 

We wish to find sub-pixel offsets, Aa; and Ay, that create 
the best match of a PSF model P to a difference image D. 
Expand the PSF model to first order in Ax and Ay, 

P = Pa + Ax + AyPy, (Al) 

and thereby write the residual imag^l as 

R = D AFP, (A2) 

where Po is the unshifted PSF, AF is the difference flux, 
and the x and y gradients of the unshifted PSF are 



dPo 

dx 



dPs, 

ay 



The PSF is normalised to 



(A3) 



(A4) 



when summed over pixels i. With = Var [Di\ the variance 
of the data Di in pixel i, the statistic measuring the 
"badness-of-fit" is 



X = 



\R[ 



Here we adopt the convenient notation 
{A\B)^y^^ \A\'^{A\A) 



(A5) 



(A6) 



for the inverse-variance weighted "dot product" of "image 
vectors" A and B, and note that x is the squared norm of 
the residual image R. 

Starting with AF = Ax = Ay = 0, we calculate the dif- 
ference flux AF, for fixed Aa; and Ay, by optimally scaling 
the shifted PSF P to fit the difference image D, giving 



AF : 



{D\P) 



and the corresponding variance 
Var [AF] = . 



(A7) 



(A8) 



Next we update the image position, for fixed AF, by scaling 
the PSF gradient images to fit the residuals. 



Ax ^ Ax + 



{R\P. 
AF\P^ 



Ay 



Ay+i^£4 



(A9) 
(AlO) 



Note that Aa; > shifts the PSF peak to smaller x, and sim- 
ilarly for y. We adopt this sign convention to simplify the equa- 
tions. 
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with variances 
Var [Ax] = 



Var [Ay] 



{AFf~\K^ 



{afYJpJ 



(All) 
(A12) 



The above results minimise = if ^F, Ax and 
Ay are independent, and if ai are fixed. As these assump- 
tions are only approximately true, iteration is required. We 
find that the iteration is faster and more stable if we take 
account of R being a linear function of Aa; and Ay. We then 
have two coupled equations, 

{D-AF (Po + AyP,)|P.) 

AF\P^i' ' 

(D-AF {Po-AxP,)[Py) 

aflpTP ' 



Write these in matrix form as 
AFAx 



H ■ 



AF Ay 
with the Hessian matrix 



(D- AFPo|P.> 

{D-AFPo[Py} 



The solution is 



\Py) 



AFAx 
AF Ay 



(£>-AFPo|Px> 
(D- AFPolP,,) 



where the inverse of the Hessian matrix is 



1 



(P. \Py) 



det{H) \ - (P, [Py) 
with the Hessian determinant 



det{H) = |P, 



- {P^ \PyY 



The sub-pixel shift is then 

{D - AFP0\{[Py[^P^ - (P. [Py) Py)) 

AF det{H) 

{D - AF P0\{\P^fPy ~ {P:,[Py) P^)) 



Ax 



Ay = 



(A13) 
(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 
(A21) 



AF det{H) 

Since is the parameter covariance matrix, the dl 
agonal elements give the variances 

D |2 

Var [Aa; 



(AF)^ det{H) 
^-'^^]- (AFfi[t(H) ' 

and the off-diagonal element gives the covariance 



Gov [Ax, Ay] 



-{P^\Py) 
{AFf det{H) 



(A22) 
(A23) 

(A24) 



Note that with AF in the denominator, these expres- 
sions become problematic when AF ~ 0. Such images carry 
very little information about the target location. Fortu- 
nately, when we optimally average over several images, the 
inverse-variance weights shift the AF factors to the numer- 
ator, so that these images receive low weight. 



A2 lens position from many images 

In fitting a microlensing dataset, we have many difference 
images Dj, and the corresponding PSFs Pj. The above 
analysis provides estimates (with error bars) that we cor- 
respondingly label AFj for the difference fluxes, Aa;j and 
Ayj for the sub-pixel offsets. 

The difference fluxes are different for each image, but 
the sub-pixel shift establishing the lens position on the ref- 
erence image should be the same for all images. The optimal 
average of the estimates Axj from individual images is 



Ax 



^ wj Aa 
i 



Var [Ax] 



(A25) 



with inverse- variance weights Wj = l/Var[Axj]. Explicit 
evaluation using (|A20[I and (|A22|) gives 



^AF (d-AFPo[P.- 



(P^ \Py) p 



Var [Ax] = 



^(AF)^ 



(P. |P 



(A26) 



(A27) 



The corresponding expressions for Ay and Var [Ay] are 
found by reversing x and y. For clarity we omit the index j 
that labels every term in the sums over images j. 

Note that the AF factors appear in the numerator only, 
so that difference images with AF^ ~ are included in the 
sums but with appropriately low weight. In our implemen- 
tation, the algorithm typically converges to ~ 10"'^ pixels 
in ~ 4 iterations. 



APPENDIX B: COMPUTING THE POISSON 
REFERENCE FLUX 

We retain here the convention (Eqn. [TJ where a difference 
image is on the same effective exposure scale as the reference 
image and a negative difference flux results when the target 
star is brighter than it is on the reference. The expected 
pixel-integrated star flux on a target image is then 



(F) = (Fo - AF) s 



(Bl) 



where Fo is the pixel-integrated flux of the lens star in the 
reference image (ADU), AF is the pixel- integrated differen- 
tial flux of the lens star in the difference image (ADU) and s 
is the exposure scale factor between the reference image and 
the target image (Eqn. [Sjl. Assuming a noiseless reference 
image, the variance in flux of the lens on a single difference 
image, is given approximately by 



Var[F] = ^^^2i^ + .^P"'^='^^ 



Fo - AF 



(B2) 



where ag is the readout noise variance (e~/pix)'^, g is the 
gain (e~/ADU), Fg^y is the background flux (ADU/pixel) 
on the target image and A^pix is the effective number of pixels 
in the photometric aperture. The different denominators for 
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each of the three terms in Eqn. IB2I are due to the read-out- 
noise being measured in units of electrons, _Fsky in units of 
target-frame ADU and Fo — AF in units of reference-frame 
ADU. An estimate of the reference flux from this single dif- 
ference image is therefore 

Fo=gs Var [F] ~ ^p'" - ^'^'^^ + Af (B3) 
gs s 

with variance 

Var [Fo] ~ Var [AF] . (B4) 

An optimal estimate for Fo is obtained by combining such 
measurements from all target images, 

j VarlFo A 

< Fo >= ^ 1 " , (B5) 
with associated variance, 

Var [< Fo >] = ^ , . (B6) 

Var[Fa,j] 

For aperture photometry measurements, A^pix is the number 
of pixels in the photometric aperture. In the case of optimal 
PSF-fitting photometry, A''pix is the effective number of sky 
pixels, which we define to be equal to the variance in the flux 
measurement that is due to the background divided by the 
background variance per pixel. To estimate A''pix, consider 
the background-limited case, where the noise variance is the 
same on each pixel, Vo, and is the dominant contributor to 
the variance in the flux measurement. If our optimal extrac- 
tion is confined to some aperture, then the variance in the 
measured flux is 



Var [AF] = ^ _ ^^^ J ^ ApixFo , (B7) 



x,y Vq 

where P{x,y) is the PSF and hence 



2 



Under a Gaussian PSF, 

Pir) = ^e-^/-^ , (B9) 
where A is the full- width at half-maximum, 

equivalent to an aperture of radius 2cr. For a Gaussian trun- 
cated at r = i?, this equates to 

»r ,_,2 (1 - e ) ttA (1 ^ ^ /-Rn-> 



